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SUMMARY 


It is shown that the neglect of the terms in the differential 
equations of atmospheric oscillations which are quadratic in the veloci- 
ties is .not justified in the case of free oscillations for heights 
greater than about 100 kilometers. This result is illustrated for two 
atmospheres of different temperature distributions where it is shown 
that at about 130 kilometers the neglected quadratic terms become equal 
to the retained linear terms and that they even exceed the latter at 
heights greater than 130 kilometers . A method is outlined for including 
the quadratic terms in the theory of atmospheric oscillations. The 
bearing of this result on atmospheric tides is discussed. 


INTRODUCTION 


In the theory of tides* whether in the ocean or in the atmosphere, 
it has been customary since the days of Laplace to neglect the terms in 
the differential equations which are quadratic in the velocities. ' This 
procedure is justified in the case of the classical theory of oceanic . 
tides, since the tidal velocities are of the same order of magnitude 
throughout the depth of the ocean and their squares can be neglected. 

In the case of the tides in the atmosphere which are similar in struc- 
ture to long waves propagated horizontally with vertical wave fronts, 
the particle velocities increase generally with height because of the 
decrease of density with height. At a sufficiently high level the 
particle velocity approaches the local velocity of sound, so that at 
such heights it is no longer permissible to neglect the squares of the 
velocities in comparison with their first powers. In the case of the - 
atmospheric oscillation' of largest amplitude, namely, the solar semi- 
diurnal tide, the particle velocity near the ground is about 35 centi- 
meters per second, which is only 0.1 percent of the velocity of sound 
there. At 100 kilometers, however, the amplitude of the particle 
velocity, as deduced from the linearized theory, is greater than at the 
ground by a factor of 100 or more, while the local sound velocity remains 
of the same order of magnitude. As long as the quadratic terms, as 



2 


NACA TN 2311; 


computed from the linearized theory, turn out to be negligible in com- 
parison with the linear terms, the linearization can be justified a 
posteriori. When, however, the quadratic terms turn out to be comparable 
with the linear terms, the linearized theory no longer applies. 

This work was conducted under the sponsorship and with the financial 
assistance of the National' Advisory Committee for Aeronautics. 


SYMBOLS 
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H, 
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u, w 


V 


X, z 

y 



acceleration due to gravity 
depth in a uniform ocean 
height of homogeneous atmosphere 

wave number 
pressure 

■unperturbed pressure at any level 

amplitude of pressure variation caused by wave 

gas constant 

temperature at any level 
time 

components of velocity 

phase velocity of wave (a/k) 

distances in x- and z-directions 

ratio of specific heats 

_ _ 2 _ ' 

" Wo 


p 


density of air 
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Pq unperturbed density 

ct frequency of wave 

ft gravitational potential of tide-producing body 

ORDER OF MAGNITUDE OF QUADRATIC TERMS 


In this section the order of magnitude of the quadratic terms, which 
are neglected in the linearized theory, will be determined. Consider a 
long wave propagated in the x-direction with a vertical wave front, and 

assume a time factor where a denotes the frequency and k, 

the wave number. For the purpose of this discussion the earth's rotation 
may be neglected and the differential equation of the principal component 
of velocity u may be written as 





+ w 


3u\ 
3z ) 


dp _ 3ft 
dx ^ dx . 


( 1 ) 


where w denotes the vertical component of velocity, p, the density, 
p, the pressure, and ft, the gravitational potential of the tide- 
producing body. In the linearized tidal theory one substitutes p Q 3u/3t 

for the term p du/dt, where p^ denotes the unperturbed density. By 


this substitution equation (l) becomes linear in u. In order to estimate 
the order of magnitude of the error committed thereby, it will suffice 
to determine whether the term u 3u/3x is negligible in comparison with 
3u/3t .Now - - 
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ft + u £ = iCTU ^ ' (k/a)u ] 
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= icru[_l - (u/v)] 



k 


NACA TN 231k 


where V denotes the phase velocity of the wave (i.e., a/k) . The ratio 
of the neglected quadratic terms in the differential equations to the 
retained linear terms is therefore of the order of the ratio of the 
particle velocity to the phase velocity of the wave. The phase velocity 
is comparable with the mean sound velocity in the vertical column of the 
atmosphere in which most of the wave energy is concentrated. Now it 
follows from the linearized equation (1) that 

Pi = p - p 0 = p 0 uV (3) 

u Pj V z)liT o( z > 

7 p 0 tf2 P 0 ( z )' f2 


with 


P 1 (z)H 0 (z) 

p q( z ) H 


(U) 




(5) 


.The factor Hq(z)/H varies relatively little with height, so that u/V 
is essentially determined by the variation of p^z^p^z), which is the 

ratio of the amplitude of the pressure variation caused by the wave to 
the unperturbed pressure at any level. 


The magnitude of u/V in equation (U) can be exhibited most easily 
for the case of the free oscillation of an atmosphere of constant 
temperature .Here , 


P 0 (z) = P 0 (0) exp (-z/H 0 ) 
P 0 ( z ) = P Q (0) exp (-z/H 0 ) 


(6) 


H = yH 0 


> 


P 1 ~ P i^ 0 ^ exp t-^ z + " kx)] 


( 7 ) 
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X = 


y - 1 
yHo 



u 

IV 


pf(o) 

yp 0 (°) 


exp 


(9z/7H 0 ) 


( 8 ) 


exp (z/5 - 6) (9) 

I 

where z is expressed in kilometers, and the values H = 7.87 kilo- 
meters and p-^(O) = 1 millimeter have been assumed to correspond to 

resonance with the solar semidiurnal tide. The right-hand side of 
equation (9) approaches unity at I 4 .O kilometers, so that above J 4 O kilo- 
meters the neglected terms in the differential equations of tidal motion 
are larger than the linear terms, which have been retained. • 


Before proceeding to a discussion of this result, the results will 
be given here of similar determinations of |u/v| for two atmospheres 
in which the temperature is not assumed to be constant with height but 
to have a distribution more closely in agreement with observations. 

for an atmosphere having the temper- 


Figure 1 gives the ratio 


u i 
dxj 6t 


the atmosphere B treated 
possess a resonance period 


ature distribution shown by curve T. This is 
by Pekeris in reference 1, which was found to 

close to 12 solar hours, with H = 8.19 kilometers. It is seen that 
the quadratic terms are less than 10 percent of the linear terms in the 
first 80 kilometers of height. At 100 kilometers this ratio reaches a 
value of 28 percent, to become greater than unity above about 123 kilo- 
meters. Figure 2 gives similar results for an atmosphere treated by ' 
Weekes and Wilkes (reference 2), which was also found to resonate with 
a period close to 12 solar hours. Here, the quadratic terms are less 
than 11 percent of the linear terms in the first 100 kilometers of height 
but exceed the linear terms at elevations greater than about 130 kilometers 


DISCUSSION OF RESULTS 


In the previous section it was shown that the linearized theory of 
atmospheric oscillations ceases to be valid at elevations greater .than 
80 to 100 kilometers, because at those heights it is no longer permis- 
sible to neglect in the differential equations the terms which are quad- 
ratic in the velocities. The inclusion of these quadratic terms in the 
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theory is a task of much greater complexity than the solution of the 
original linearized equations. A procedure which could be followed in 
obtaining a solution of the complete differential equations is to expand 
the exact solution into a series proceeding in powers of the amplitude, 
starting with the solution of the linearized equations as a first approxi- 
mation. Such an expansion would be expected to converge rapidly below 
about 100 kilometers. At higher elevations the exact solution would be 
expected to differ materially from that represented by the first 

approximation. ^ 

Pending the development of the exact solution for atmospheres of 
given temperature distributions, a brief qualitative discussion will now 
be given of the possible effects of the neglected quadratic terms in the 
tidal equations on the resonance theory of atmospheric tides. 

Higher harmonics .- The solution of the exact wave equations for the 

free oscillations of the atmosphere will be periodic in time with an 
admixture of higher harmonics whose relative magnitude increases with 
elevation. A sinusoidal tidal potential will excite in the lower atmos- 
phere a nearly pure sinusoidal oscillation of the same period, but at 
high elevations the oscillation will contain, besides the fundamental, 
also strong higher harmonics. The amplitude of the fundamental in the 
E layer may differ considerably from that predicted by the linearized 
solution, because the quadratic terms in the equations also contribute 
to the fundamental. 

Free period.- The quadratic terms will cause a change in the free 

period of oscillation of the atmosphere or, more directly, in the velocity 
of propagation of long waves. This change would, however, be expected to 
be small in those atmospheres in which the linearized solution shows that 
only a small fraction of the wave energy in an atmospheric column is 
contained in the upper layers where the quadratic terms become significant. 
In the example of an atmosphere of constant temperature treated in the 
preceding section, the quadratic terms become dominant above I 4 .O kilometers. 
The part of the atmosphere above J 4 .O kilometers carries, however, only 
10 percent of the total wave energy in the atmospheric column, so that it 
is not likely that the inclusion of the quadratic terms will appreciably 
alter the period computed from the linear theory. This, would also hold 
for the atmosphere (B) of figure 1, since it is seen that the area under 

^Some exploratory work on this method of solution has been made, 
showing that on account of the earth's rotation the second harmonic will 
not be amplified by resonance, as is the first harmonic. A contrary 
result would have proved a serious difficulty in the proposed method. 
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the E curve (energy density) above the 80 -kilometer level is a small 
fraction of the total area under this curve. On the other hand, if the 
temperature distribution in the atmosphere is such that a substantial* 
portion of the wave energy in a column is located at the high levels 
where the quadratic terms are significant, then the free period of such 
an atmosphere may- differ radically from that computed on the basis of 
linear theory. An atmosphere of the latter type is one' having a positive 
temperature lapse rate in the outer layers. 

Resonance theory.— It would follow from the preceding discussion 

that, provided the temperature distribution in the outer layer of the 
atmosphere is of the type shown in figure 1 , the linear theory yields a 
good approximation to the amplification of the tide and the free period. 

Thus the claims of the resonance theory are not weakened for this type of 
atmosphere by consideration of the effects of the quadratic terms. On 
the other hand, the linear theory cannot be relied upon for a description 
of the features of the tidal oscillations in and above the E layer. 

Effects of temperature distribution above 100 kilometer s.- Since in 
the E layer and higher elevations the quadratic terms in the tidal equations 
become dominant, it is clear that any consequences, drawn on the basis of 
the linear theory , which result from assumed specific changes in the tem- 
perature distribution in and above the E layer are open to question. Thus 
the explanation given by Weekes and Wilkes (reference 2) of the observa- 
tions on the lunar tide in the E layer made by Appleton and Weekes (refer- 
ence 3 ) cannot be accepted, because it is dependent on the existence of 
a temperature maximum in the E layer where the quadratic terms turn out 
to be not only not negligible but twice as large as the retained linear 
terms. 

The question as to whether the resonance theory of atmospheric tides 
is reconcilable with an assumed positive temperature lapse rate in the 
outer atmospheric layers must await a detailed attack on the solution of 
the nonlinearized tidal equations. 

Effect of density distribution on atmospheri c oscillations.- The 

emergence of the quadratic terms in the. tidal equations is but a concealed 
effect of the density on the atmospheric oscillations, since the increase 
with height of the amplitude of the particle velocities is due to the 
decrease with height of the density. 


Institute for Advanced Study 

Princeton N. J., September 22, 19£o 
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Figure !.- Ratio of neglected term u(8u/dx) to term au/at which is 
retained in equations of atmospheric oscillations. Curve T is ass um ed 
temperature distribution; E, the wave-energy density. 
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Figure 2.- Ratio of neglected term u(3u/3x) to term 3u/3t which is 

retained in theory of atmospheric oscillations. Curve T shows assumed 
temperature distribution. 
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